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ABSTRACT 

The large-eddy simulation of the spatial evolution of a station- 
ary crossflow vortex packet in a three-dimensional boundary layer 
was performed. Although a coarse grid was used (compared to that 
required by a direct numerical simulation) the essential features of 
the disturbance evolution, such as the spanwise disturbance spread- 
ing and the vortex rollover, were captured accurately. The eddy 
viscosity became significant only in the late nonlinear stages of the 
simulation. 

1. Introduction 

The understanding and control of the transition from laminar to 
turbulent flow on swept wings are of great technological importance 
for the aerospace industry for many reasons. Laminar flow control, 
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for instance, would reduce the drag and increase the efficiency of 
transportation systems. 

The transitional flow on swept wings, however, involves extremely 
complex phenomena. Leading-edge contamination and crossflow in- 
stabilities, for example, can cause the flow to breakdown. Arising at 
the wing/fuselage juncture, turbulence contaminates the flow along 
the leading edge, causing transition. Gaster (1965) showed that this 
turbulent contamination can be avoided by placing a bump on the 
leading edge and near the wing/fuselage juncture. Crossflow insta- 
bilities are excited because of the pressure gradients and crossflow 
velocity component, which develops because of wing sweep. The 
crossflow velocity has an inflectional point leading to an inviscid in- 
stability, which is very different from the viscous instability of the 
T-S type. Among others, Muller and Bippes (1988) have observed 
both stationary and traveling crossflow disturbances in their wind- 
tunnel experiments. 

In experimental studies of transition, however, such factors as 
free-stream turbulence, wind-tunnel noise, wall interference, surface 
roughness, measurement limitations, etc. make the identification of 
the instabilities that are of primary importance a very difficult task. 
Theories also have difficulties when the nonlinear interaction of dis- 
turbance is important. Numerical simulations which solve the un- 
steady, nonlinear Navier- Stokes equations provide a very promising 
but expensive way to investigate the three-dimensional boundary- 
layer transition problem. The direct numerical simulations (DNS) 
to date have been limited to very low Reynolds numbers because 
of the grid-resolution requirements; for example, Spalart (1989) per- 
formed a DNS of swept Hiemenz flow with a 256 X 20 X 32 mesh 
(chordwisex wall-normal X span wise). The evolution of the distur- 
bance up to the very weakly nonlinear stage was captured in the 
simulation. With DNS, Joslin and Street (1993) studied the role 
of stationary crossflow vortices on swept-wedge flow; using a grid 
of 901 X 61 X 32 points, they captured the qualitative breakdown 
features of nonlinear stationary vortex interactions. 

A technique that could quantify the fuU transition process on a 
swept wing is the large-eddy simulation (LES) of the Navier-Stokes 
equations. In LES the large scales of motion are computed accu- 
rately, while only the subgrid- scale motions, which are presumably 
more universal, are modeled. This permits the use of simpler parame- 
terizations of the subgrid-scale (SGS) stresses. Recently, for instance, 



the dynamic eddy- viscosity model has been successfully applied to 
the plane channel transition problem by Germano et al . , 1991. The 
most attractive characteristic of the model is that the model coef- 
ficient is calculated from the large-scale flow features rather than 
specified a priori ; this feature allows the SGS stresses to vanish in 
laminar flows and near solid walls. 

In this paper, LES is used to study the role of stationary crossflow 
vortices in swept-wing boundary-layer flow. Section 2 gives the the 
problem formulation, section 3 a discussion of the results, and section 
4 presents some conclusions. 


2. Problem Formulation 

The present technique relies on the decomposition of the instan- 
taneous velocity fq(x,t) and pressure p(x,/) into base, P;(x) and 
P(x), and disturbance, rq(x,t) and p(x,/), components as 

Ui(x.,t) = Ui(x.) + Ui(x.,t) and p(x, t) = P(x) + p(x, t), (1) 


where x = ( x,y,z ) are chordwise, wall-normal, and spanwise coordi- 
nates and t is time. 

After substituting (1) into the Navier- Stokes equations and sub- 
tracting out the base-flow equations, the disturbance equations result 
and are given as 
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In LES, the large-scale (grid-resolved) components of the velocity 
and pressure are calculated and the effects of the small, unresolved 
scales are modeled. By applying the filtering operation 
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to (2) and (3), the governing equations for the large-scale velocity 
and pressure can be obtained: 
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and 
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where ( ) denotes the filtered variables (corresponding to the large- 
scale motions) and is the subgrid-scale (SGS) stress tensor given 

by 

Tij = V^Uj -UiUj, (7) 

Here, the tensor is modeled with the dynamic SGS viscosity model 
proposed by Germano et al. (1991) in the later formulation described 
by Lilly (1992). 

To solve equations (5)-(7), fourth-order finite differences and fourth- 
order compact differences are used in the chordwise direction for the 
pressure and the momentum equations, respectively, Chebyshev se- 
ries are used in the wall-normal direction, and Fourier series are 
used in the spanwise direction. The implicit Crank- Nicolson scheme 
is used for the time- advancement of the wall-normal diffusion terms 
and a three-stage Runge-Kutta scheme for the remaining terms. 

A schematic of the swept-wing problem is shown in Figure 1. 
To simplify the problem, wing curvature is neglected and an infinite 
span is assumed. For the resulting swept-wedge problem, the base 
flow can be described by the Falkner-Skan- Cooke similarity profiles. 
Furthermore, as shown in Figure 1, the computational domain is only 
a chordwise subset of the airfoil. 

The disturbances are required to vanish far from the wing and 
on the wing surface, the steady base flow is forced at the inflow 
boundary, and the buffer-domain technique of Streett and Macaraeg 
(1989) is used at the outflow boundary. Periodic boundary conditions 
are implied in the spanwise direction because Fourier series are used. 
Stationary crossflow- vortex packets are introduced into the boundary 
layer by steady suction and blowing through a spanwise periodic 
array of holes in the wing. The suction and blowing distribution is 
shown in Figure 2. 

The parameters for the simulation are chosen to match those of 
wind-tunnel experiments (Muller and Bippes 1988; Dagenhart et al. 
1989). The chord length is c = 1.84m, the wing sweep is q ^ = 45°, 
and the angle- of- at tack is approximately —4°. Based on the inflow 
boundary-layer thickness and edge velocity, the Reynolds number is 
Re = 996.4. The disturbance is initiated near x c = 0.22, where x c is 
the relative chordwise coordinate x/c. 




Figure 1: Schematic of the swept wing, computational domain and 
coordinate directions. 
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Figure 2: Suction/blowing velocity profiles in the chordwise and 
spanwise directions 


The computational domain is 8578 X 50tf X 108tf (respectively in 
the chordwise, wall-normal and spanwise directions), and 301x41x33 
grid points are used for the LES. The previous DNS by Joslin and 
Streett (1993) used 901 X 61 X 33 grid points in the same domain; 
they note that the resolution appeared marginal in the late nonlinear 
stages of breakdown. Time was advanced with a fixed step size At = 
0.4, and the amplitude of the suction and blowing velocity profile 
(Figure 2) was A = 1 X 10 -5 . 

3. Results 

Figure 3 shows instantaneous (base+disturbance) velocity pro- 
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Figure 3: Total (base+disturbance) velocity profiles at different 

cliordwise and span wise locations. U; U + u at z =: o 0> 

□ 3.375, A 6.75, o 10.125, * 13.5, X 16.875, + 20.25, • 23.625. 
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Figure 4: Cliordwise disturbance velocity (u) contours at the chord- 
wise locations: (a) ,r c =0.25, (b) ,r c =0.35, and (c) ,r c =0.45. 

hies obtained from the LES at two stations along the chord. The 
various profiles at each station correspond to adjacent spanwise lo- 
cations. The cliordwise and spanwise velocity profiles indicate that, 
at x c = 0.375, the growing disturbance has already altered the mean 
how. Downstream, at x c = 0.425, the amplitude of the disturbance 
velocities become sufficiently large to cause inflections to the velocity 
profiles. These inflectional profiles were observed previously in both 
DNS (Joslin et al ., 1993) and experiments (Muller and Bippes 1988, 
Dagenhart et al. 1989). The normal velocity profiles show evidence 
of spanwise regions of how moving away from the wall and adjacent 
spanwise regions of how moving toward the wall; these types of mo- 
tion intensify as one moves downstream. Although not shown in this 
study, these inflectional profiles spawn the rapid growth of secondary 
instability modes, which lead to the catastrophic breakdown to tur- 
bulence (see Koliama et al ., 1991). The quantitative agreement with 
the DNS results of Joslin and Streett (1993) is quite good. 

Figure 4 shows spanwise planes of the cliordwise disturbance ve- 
locity contours obtained by the LES at three cliordwise locations. 
At x c = 0.25 (near the suction and blowing hole), the vortex packet 







evolves isolated from adjacent packets. Downstream near x c = 0.35, 
the disturbance spreads in the spanwise direction and interacts with 
adjacent packets. Finally, near x c = 0.45, the crossflow vortices non- 
linearly interact causing the so-called vortex- rollover phenomenon. 
This phenomenon is characterized by low-speed fluid near the wall 
being lifted up into the boundary layer while high-speed fluid is 
drawn toward the wall. During the stage of vortex- rollover the in- 
flectional velocity profiles discussed above are observed. Once again, 
the LES are in good agreement with the DNS. 

Figure 5 shows the maximum amplitude of the disturbance as a 
function of the cliordwise position. The disturbance is found to grow 
exponentially, with a sudden increase in the growth rate near the 
spanwise spreading and vortex/ vortex interaction region (x c = 0.32). 





Figure 6 shows that the maximum ratio of eddy to molecular viscos- 
ity remains very small until the nonlinear vortex/ vortex interactions 
begin. As the disturbance amplitudes grow, the magnitude of the 
eddy viscosity grows and the SGS stresses begin to become impor- 
tant. Note that no ad hoc adjustment was made to switch the eddy 
viscosity on or off, and its increase in magnitude reflects the capa- 
bility of the dynamic model to adjust to the characteristics of the 
disturbed flow held. 


4. Conclusions 

The dynamic subgrid scale eddy- viscosity model was used for the 
large-eddy simulation of swept-wing transition. The eddy viscosity 
remained small throughout most of the computational domain, but 
started to grow in the region of strong nonlinear vortex/vortex in- 
teractions, adjusting to the local state of the how. Conceivably, if 
unsteady modes were present in the how held and the simulation were 
continued, the role of the eddy viscosity in damping the unresolved 
small-scale huctuations would become dominant and the inhectional 
prohles would cause secondary instabilities to grow at an explosive 
rate, causing transition. 

Because the magnitude of the eddy viscosity was small in most of 
the computational domain, essentially a coarse-grid DNS captured 
the essential features of the spatial evolution of the crosshow vortex 
packet. Signihcant among these features are the vortex- rollover and 
the subsequent generation of inhectional velocity prohles, which are 
responsible for the rapid transition to turbulent how. The present 
results were in good agreement with DNS data. 

Based on the preliminary results, it appears that LES can be a 
useful tool to quantify the phenomenon in the later stages of tran- 
sition on a swept wing, where the resolution required by a DNS be- 
comes prohibitively expensive. Efforts along this route are currently 
underway. 
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